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Some of the parameters we call "constants of Nature" may in fact be 
variables related to the local values of some dynamical fields. During infla- 
tion, these variables are randomized by quantum fluctuations. In cases when 
' the variable in question (call it x) takes values in a continuous range, all 

thermalized regions in the universe are statistically equivalent, and a gauge 
^ - invariant procedure for calculating the probability distribution for x is known. 

. This is the so-called "spherical cutoff method" . In order to find the probability 

distribution for x it suffices to consider a large spherical patch in a single ther- 
malized region. Here, we generalize this method to the case when the range of 
X is discontinuous and there are several different types of thermalized region. 
We first formulate a set of requirements that any such generalization should 
satisfy, and then introduce a prescription that meets all the requirements. We 
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bubble universes in the open inflation scenario. 



I. INTRODUCTION 



The parameters we call "constants of Nature" may in fact be variables related to the 
local values of certain dynamical fields. For example, what we perceive as a cosmological 
constant could be a potential U{x) of some slowly varying field x{^)- If this potential is very 
flat, so that the evolution of x is much slower than the Hubble expansion, then observations 
will not distinguish between U{x) ^-iid a true cosmological constant. Observers in different 
5^ \ parts of the universe could then measure different values of Uix)- 

Spatial variation of the fields Xa associated with the "constants" can naturally arise 
in the framework of inflationary cosmology [|1|. The dynamics of light scalar fields during 
inflation are strongly influenced by quantum fluctuations, so different regions of the universe 
thermalize with different values of Xa- An important question is whether or not we can 
predict the values of the "constants" we are most likely to observe. In more general terms, 
we are interested in determining the probability distribution 'P(x) for us to measure certain 
values of Xa- The answer to this question must involve anthropic considerations to some 
extent. The laws of physics may be sufficient to determine the range and even the spacetime 
distribution of the variables Xa- However, some values of Xa which are physically allowed 
may be incompatible with the very existence of observers, and in this case they will never 
be measured. The relevant question is then how to assign a weight to this selection effect. 

The inflationary scenario implies a very large universe inhabited by numerous civilizations 
that will measure different values of Xa- We can define the probability V{x)dxi---dxk for 
Xa to be in the intervals dxa as being proportional to the number of civilizations which will 



measure Xa in that interval 0. This includes all present, past and future civilizations; in 
other words, it is the number of civilizations throughout the entire spacetime, rather than 
at a particular moment of time. Assuming that we are a typical civilization, we can expect 
to observe Xa near the maximum of 'P(x) 0- The assumption of being typical has been 
called the "principle of mediocrity" in Ref. 0. 

An immediate objection to this approach is that we are ignorant about the origin of life, 
let alone intelligence, and therefore the number of civilizations cannot be calculated. But 
even if this were true, the approach can still be used to find the probability distribution for 
parameters which do not affect the physical processes involved in the evolution of life. The 
cosmological constant A, the density parameter Q and the amplitude of density fluctuations 
Q are examples of such parameters. Assuming that our fields Xa belong to this category, the 
probability for a civilization to evolve on a suitable planet is then independent of Xa, and 
instead of the number of civilizations we can use the number of habitable planets or, as a 
rough approximation, the number of galaxies. Thus, we can write 



where dM' is the number of galaxies that are going to be formed in regions where Xa take 
values in the intervals dxa- 

The probability distribution (|l]) based on plain galaxy counting is interesting in its own 
right, since it gives a quantitative characterization of the large scale properties of the uni- 
verse. Thus, the general rules for calculating (|l|) are worth investigating quite independently 
from anthropic considerations. These considerations can always be included a posteriori, as 
an additional factor giving the number of civilizations per galaxy. 

The number of galaxies dM{x) in Eq. (|1]) is proportional to the volume of the comoving 
regions where Xa take specified values and to the density of galaxies in those regions. The 
volumes and the densities can be evaluated at any time, as long as we include both galaxies 
that were formed in the past and those that are going to be formed in the future. It is 
convenient to evaluate the volumes and the densities at the time when inflation ends and 
vacuum energy thermalizes, that is, on the thermalization surface S*. Then we can write 



Here, V*{x)d^X is proportional to the volume of thermalized regions where Xa take values 
in the intervals dxa-, and z/(x) is the number of galaxies that form per unit thermalized 
volume with cosmological parameters specified by the values of Xa- The calculation of z/(x) 
is a standard astrophysical problem which is completely unrelated to the calculation of the 
volume factor V^:{x)i and which does not pose difficulties of principle. 
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FIG. 1. Symmetric double well inflaton potential 

The meaning of Eq. (0) is unambiguous in models where the total number of galaxies 
in the universe is finite. Otherwise, one has to introduce some cutoff and define the ra- 
tio of probabilities for the intervals dJ^x^^^ ^^^^ dP'x^'^^ s-s the ratio of the galaxy numbers 
dJ\f^^^ / dj\f^'^^ in the limit when the cutoff is removed. However, this limiting procedure has 
proved to be rather non-trivial, and a general method which would apply to all possible 
eternally inflating scenarios has not yet been found. 

The situation is relatively straightforward in the case of an infinite universe which is 
more or less homogeneous on very large scales. One can evaluate the ratio dj\f^^y dj\f^'^^ in 
a large comoving volume V and then take the limit as V — > oo. The result is expected to be 
independent of the limiting procedure; for example, it should not depend on the shape of 
the volume V. (It is assumed that the volume selection is unbiased, that is, that the volume 
V is not carved to favor some values of Xa at the expense of other values.) 

However, the situation with an infinite universe which is homogeneous on very large 
scales is not generic in the context of inflation. Most inflationary scenarios predict that 
inflation is eternal to the future, and therefore the universe is never completely thermalized 
(for a recent review of eternal inflation see An example which is particularly 

relevant to the subject of the present paper is given by the double well inflaton potential 
depicted in Fig. 1. The inflaton can thermalize in two different vacua, labeled by rji and 
ri2. The spacetime distribution of the field in this model is depicted in Fig. 2. There are 
thermalized regions of two types, characterized by the inflaton vacuum expectation value r)i. 
Thermalized regions with (p = rji are disconnected from thermalized regions with (j) = r]2- 



Both types of regions are separated by inflating domain walls |13,14], and so the universe is 



never completely thermalized. Each thermalization surface is infinite, so it will contain an 
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infinite number of galaxies. Moreover, tliere are an infinite number of tliermalized regions 
of eacli type. Tlierefore, tlie implementation of Eq. (|l]) for calculating probabilities requires 
the comparison of infinite sets of galaxies which lie in disconnected regions of the universe. 
A similar spacetime structure is obtained if we make the inflaton potential periodic by 
identifying the two minima. In this case there are still two different types of thermalization 
surfaces, characterized by the two topologically different paths that one can take from the 
top of the potential to the thermalized region. Although the particle physics parameters 
of both types of thermalized regions are guaranteed to be the same, other cosmological 
parameters such as the spectrum of density perturbations or the spatial distribution of an 
effective cosmological constant will in general be different. 




FIG. 2. A numerical simulation of the spacetime structure of an inflating universe [11|. The 
simulation corresponds to a double-well inflaton potential, with two degenerate minima where the 
inflaton takes the values iry. Inflating regions are white, while thermalized regions with inflaton 
values equal to +ry and —r] are shown with different shades of grey . 



In any model of eternal inflation, the volumes of both inflating and thermalized regions 
grow exponentially with time and the number of galaxies grows without bound, even in a 
region of a finite comoving size. One can try to deal with this problem by introducing a time 
cutoff and including only regions that thermalized prior to some moment of time tc, with 
the limit tc ^ oo at the end. One finds, however, that the resulting probability distributions 
are extremely sensitive to the choice of the time coordinate t 0]. Coordinates in General 
Relativity are arbitrary labels, and such gauge dependence of the results casts doubt on any 
conclusion reached using this approach. 

A resolution of the gauge dependence problem was proposed in Ref. [110] and subsequently 
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developed in [|TT|. The proposed method can be summarized as follows. Let us first assume 
that infiating and thermalized regions of spacetime are separated by a single thermalization 
surface S^,. The problem with the constant-time cutoff procedures is that they cut the 
surface T,^ in a biased way, favoring certain values of x ^i-nd disfavoring other values. We 
thus need a portion of T,^ selected without bias. The simplest strategy is to use a "spherical" 
cutoff. Choose an arbitrary point P on S*. Define a sphere of radius R to include all points 
Q whose distance from P along is d{Q, P) < R. We can use Eq. ([^) to evaluate the 
probability distribution 'P(x) in a spherical volume of radius Rc and then let Rc oo- If 
the fields Xa vary in a finite range, they will run through all of their values many times in 
a spherical volume of sufficiently large radius. We expect, therefore, that the distribution 
V{x) will rapidly converge as the cutoff radius Rc is increased. We expect also that the 
resulting distribution will be independent of the choice of point P which serves as the center 
of the sphere. The same procedure can be used for fields with an infinite range of variation, 
provided that the probability distributions for Xa are concentrated within a finite range, 
with a negligible probability of finding Xa very far away from that range. 

Suppose now that there is an infinite number of disconnected thermalization surfaces, 
as it happens generically in eternal infiation. Further, we assume that the variables Xa of 
interest are such that their whole range of values is allowed to occur in a single thermalized 
region (this is the case, for instance, for the slowly varying field x which plays the role of 
a cosmological constant), and that, unlike the case of the double well potential in Fig. 1, 
there is only one type of thermalized region. We can then pick an arbitrary connected 
component of T,^, and apply the spherical cutoff prescription described above. Since the 
inflationary dynamics of the flelds Xa have a stochastic nature, the distributions of Xa on 
different connected components of S^, should be statistically equivalent, and the resulting 
probability distribution V{x) should be the same for all components. This has been verifled 
both analytically and numerically in ||11|| . 

The main shortcoming of the spherical cutoff prescription is that as it stands it cannot 
be applied to models where the inflaton potential has a discrete set of minima, as in the 
example shown in Fig. 1. More precisely, the problem arises when the minima are separated 
by inflating domain walls |T3|,|1^. In this case, we can introduce a discrete variable n 
labeling different minima. Each connected component of the thermalization surface S^, 
will be characterized by a single value of n (unless the different minima can be separated by 
non-inflating domain walls) and it is clear that the probability distribution for n cannot be 
determined by studying one such component. 

The purpose of this paper is to propose a generahzation of the spherical cut-off prescrip- 
tion that would be applicable in the general case. We begin in Section II by formulating the 
requirements that we believe any such proposal should satisfy. We require that it should be 
gauge-independent and should reduce to the spherical cutoff prescription in the absence of 
discrete variables. Moreover, we consider a class of asymmetric double-well potentials for 
which the probabilities can be calculated in a well-motivated way. We then require that the 
general prescription should give the same result for this class of potentials. In Section III we 
propose a prescription that satisfles all of the above requirements. We use this prescription 
in Section IV to calculate probabilities for bubble universes in open inflation scenario. Our 
conclusions are summarised and discussed in Section V. 
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II. REQUIREMENTS 



Suppose that the inflaton potential has minima labeled by z = 1, 2, A^, so that there 
are different types of thermalized regions. Suppose also that there is a set of scalar fields 
Xa which take a continuous range of values in all types of thermalized regions. Our goal is 
to calculate the probability distribution 

v^{x) = PAix)- (3) 

Here, 'Pj(x) is the normalized distribution for x within n-th type of thermalized region, 

/ d'xnx) = 1, (4) 

and Pi is the probability for an observer to find herself in a thermalized region of type i. 

We begin with the obvious requirements that the probabilities (§) should be gauge- 
invariant and should satisfy the normalization condition 

EP, = 1. (5) 

i 

Next, we require that for = 1 the prescription should reduce to the spherical cutoff 
prescription. This implies that the spherical cutoff should be used for the calculation of the 
individual distributions Vi{x) within each type of thermalized region. What remains to be 
determined are the relative probabilities of different types of regions, Pj. 




FIG. 3. Asymmetric double well potential which is symmetric near the maximum in the range 

-<f>s < <p < (t)s- 
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Finally, we introduce a class of inflaton potentials for which we believe there is a well- 
motivated answer for the probabilities. Consider first a symmetric double-well potential, 
V{(f)) = V{—(j)), with a maximum at = and two minima at 0i 2 = ±?7, as in Fig. 1. 
Clearly, in this case the symmetry of the problem dictates that Pi = P2 = 0.5. 

Next, we consider an asymmetric double-well, which however is symmetric in some range 
of (p near the maximum, < (ps (see Fig. 3). Quantum fluctuations of (p dominate the 
dynamics in the range |0| ^ (pq, while for (p ^ (pq the evolution is essentially deterministic. 
We shall assume that (pg is in the deterministic slow-roll regime, (pg ^ (pq. Let us consider 
constant-0 surfaces (p = ±0^. These are infinite spacelike surfaces which have \(p\ < (ps 
everywhere in their past. Since the potential V{(p) is symmetric in this range of 0, all 
these surfaces are statistically equivalent. The symmetry of the problem suggests that the 
probabilities Pi,2 can be calculated by sampling comoving spheres of equal radius on the two 
types of surfaces. The ratio of the probabilities will then be 

P1IP2 = N1/N2, (6) 

where A^^i , N2 are average numbers of galaxies that will form in large comoving spherical 
regions which have equal radii at = ±0^. If we assume for simplicity that the two types of 
regions have identical physics at thermalization and afterwards, then the difference between 
A^i and N2 can only be due to the different inflationary expansion factors Zi characterizing 
the evolution from ±0^ to the thermalization points 0*j. We then have Ni oc Zf and 

P1/P2 = {Zi/Z^f. (7) 

We shall require that the general prescription for probabilities should reproduce Eq. ([^ in 
the case of asymmetric double-well potentials of the type we discussed here. 



III. THE PROPOSAL 

In the general case, the inflaton potential has no symmetries to guide our selection of the 
equal-0 surfaces on which to calculate probabilities. Suppose the potential has a maximum, 
which we choose to be at = 0, and two minima with thermalization points at 0*i and 0*2- 
We can then choose some arbitrary values 0i and 02 in the slow roll ranges of adjacent 
to 0*1 and 0^,2, respectively, and calculate the probabilities by sampling the surfaces = 0^. 
Imagine for a moment that the number of thermalized regions of both types and the number 
of galaxies in each region are all finite. Then we could write 

^ = (S) 

P2 P2N2 

Here, Pi is the probability for a randomly selected thermalized region to be of type i and Ni 
is the average number of galaxies in a type-i region. In our case, however, the number of 
thermalized regions and the number of galaxies in each region are infinite, so the definitions 
of the probabilities pi and of the ratio N1/N2 are problematic. 

It has been remarked [^,^ that the problem of determining pi is similar to the problem 
of calculating the probability peven that a randomly selected integer is even. If we take a 
long stretch of the natural series 
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1,2,3,..., 



(9) 



it will have nearly equal quantities of even and odd numbers, suggesting that peven = 1/2. 
However, the series can be reordered as 

1,2,4,3,6,8,5,..., (10) 

and the same calculation would give peven = 2/3. It is clear that by appropriately ordering 
the series one can obtain any answer for peven between and 1. This seems to suggest that 
the probabilities Pi are hopelessly ill-defined. 

We note, however, that the situation with the natural series is not as bad as it may 
seem. The series has a natural ordering in which the nearest neighbors of each number differ 
from that number by 1, and we can require that our sampling procedure should respect this 
natural "topology". Then we have peven = 1/2, which is the answer that one intuitively 
expects. With an infinite number of thermalized regions of different types, one could also 
order the list of regions in a way that would give any desired result for Pi/p2. But again 
one can hope that this ambiguity can be removed if we require that our sampling procedure 
should reflect the spatial distribution of the regions. 

The ratio N1/N2 can be calculated by counting galaxies in spheres located in regions of 
the two types. In the double-well example of Section II, there is complete symmetry between 
the surfaces = 0s and = — 0^. We expect, therefore, that pi = p2 and we choose the radii 
of the spheres to be equal at = ±0^. However, it is not clear what sets the relative size of 
the spheres in the absence of symmetry. We shall now describe the method we propose for 
evaluating pi and N1/N2 in the general case. 

We start by noting that there is one thing that thermalized regions of the two types have 
in common. In their past they all went through a period of stochastic inflation, with the 
infiaton field undergoing a random walk near the top of the potential. Our idea is to use 
some markers from this early era for the calculation of probabilities. 

Let us imagine that markers are point objects and that they are produced at a constant 
rate per unit spacetime volume in inflating regions where is at the top of its potential, 
V{(f)) ~ Vmax- After that, the markers evolve as comoving test particles and eventually end 
up in a thermalized region of one type or the other. We shall define Pi as the fraction of 
markers that end up in regions of type i. Furthermore, the reference length scale on which 
the number of galaxies is counted in a region of type i will be set by the average distace di 
between markers in that type of region. In other words, the galaxies are counted in spheres 
of radii -Ri and R2 such that R1/R2 = di/d2. To summarize, we propose that the relative 
probabilities for the constants of nature are given by (|^) with 

^ = ^^. (11) 

P2 P2ni V2 

where rii oc d~'^ is the mean number density of markers in region of type % and i^j, given by 

is the mean density of galaxies in that type of region. As mentioned in the introduction, the 
calculation of z/j is a standard astrophysical problem which we shall not be concerned with 
in this paper. 
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A physical counterpart to the production of our ideahzed markers is quantum nucle- 



ation of black holes in an inflating universe. Black holes are produced at a rate |T5 
r oc exp[— 1/8V^(0)] which grows exponentially with V{(f)), so that by far the highest rate 
is achieved at V{(f)) = Vmax- Here, we prefer not to identify markers with black holes and 
to think of them as of idealized point objects. This frees us from concerns such as the con- 
tribution to black hole production from regions which are not at the top of the potential, 
black hole evaporation, etc. The reality and observability of markers is not really an issue, 
since the comparison of causally disconnected thermalized regions is certainly a gedanken 
experiment. 

The average comoving distance between the markers in a thermalized region is manifestly 



gauge invariant, and so is the ratio n2/ni in Eq. (|TTD . The quantities pi are also gauge- 
independent, and thus the probabilities Pi in ( |Tl| ) should be gauge-invariant. 

It is also easy to verify that the above proposal gives the expected result (J^) for the 
probabilities when the infiaton potential is symmetric in the diffusion range of 0. In this case, 
the average distance between markers on the surface = 0^ is the same as that on = — 0^ 
(with 0s defined in Section II) . Hence, the distances between markers at thermalization may 
differ only due to the difference in the expansion factors from ±0^ to 0=Kj, 

di/d2 = Z^IZ2, (12) 

and Eq.(^ follows immediately. 




FIG. 4. Asymmetric infiaton potential with a metastable vacuum at the top, which replaces 
the quantum diffusion region of Fig. 3. 
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IV. OPEN INFLATION 



As an illustration of the method we shall now calculate probabilities for a model of "open 



inflation" [Il6|-[T8[]. We assume an infiaton potential of the form shown in Fig. 4. It has a 
metastable false vacuum at = which is separated by potential barriers from two slow roll 
regions on the left and on the right. The false vacuum decays through bubble nucleation, 
and the infiaton rolls towards the true vacuum inside the bubbles. Comoving observers 
inside each bubble would, after thermalization of the infiaton, see themselves in an open 
homogeneous universe. (Hence the name "open inflation".) If the bubble nucleation rate 
is not too high, bubble collisions are rare and inflation is eternal. Assuming that the two 
types of bubbles have identical low-energy physics, we would like to find the probability for 
an observer to find herself in one type of bubble or the other. 
In false vacuum regions outside bubbles the metric is de Sitter, 

ds' = dt' - e^^^'d^^ (13) 

where Hq = SnGVo/S and Vq = ^(0) is the false vacuum energy density. The bubble interior 
has the geometry of an open Robertson- Walker universe, 

ds^ = dr^ - a' (r) [d^^ + sinh' ^dn^] . ( 14) 

After nucleation, the bubble wall expands rapidly approaching the speed of light. If the 
initial bubble size is much smaller than the de Sitter horizon Hq^, then the wall worldsheet 
is well approximated by the future light cone of the center of spacetime symmetry of the 
bubble. We shall choose coordinates so that this center is at x = t = and r = 0. Then, as 
t ^ oo, the bubble wall asymptotically approaches |x| — Hq^. The spacetime geometry of 
the bubble is illustrated in Fig. 5. 

The relation between the coordinates (t,x) and (t, ^) can be easily found if we assume, 
as it is usually done, that (i) the potential V{(j)) immediately outside the barrier has nearly 
the same value Vq as in the false vacuum, and (ii) that the gravitational effect of the bubble 
wall is negligible. Then, at sufficiently small values of r, the metric (|1^) inside the bubble is 
close to the open de Sitter metric, a(r) = Hq^ smh.{HQT) . The coordinates (t,x) and (r, ^) 
are related by the usual transformation between the flat and open de Sitter charts, which 
we shall not reproduce here. 

We assume that markers are produced at a constant rate r in the false vacuum outside 
the bubbles. The density of markers n(x, t) satisfles the equation 

— + 3Hon = r, (15) 

which has a stationary solution hq = t/SHq. This de Sitter-invariant solution is rapidly 
approached, regardless of the initial conditions. We are interested in the average separation 
between the markers on the surfaces r = r,,j corresponding to thermalization of the inflaton 
in the two types of bubbles. The densities of markers n^i on these surfaces are also constant, 
due to the spacetime symmetry of the bubbles (all points on the surface r = const are 
equivalent). 
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FIG. 5. Spacetime geometry of the bubble in co-moving coordinates. Here, r = |x|, and various 
r = const, surfaces of the open Robertson- Walker universe inside the bubble are shown, as well as 
one of the ^ = const, curves. 

The interior geometries of different types of bubbles are nearly identical at early times 
(small r), while V{(j)) ~ Vq. Let us choose a value tq in that range and consider surfaces 
T = tq. The geometry and the distribution of markers in the past of such surfaces are 
the same for the two types of bubbles, and therefore the average separations of markers on 
these surfaces should also be the same, cii(ro) = ci2(T"o) = do. The separation of markers at 
thermalization is = Z^ido, where Z^i is the expansion factor between tq and r^,i, and we 
have 



To complete the calculation, we have to determine the fraction pi of the markers that 
end up in type-z bubbles. Let fi(t) be the fraction of comoving volume occupied by type-i 
bubbles in a comoving volume which we assume to be free of bubbles at the initial moment 




(16) 



t = 0, 




(17) 




(18) 




(19) 
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where Aj is the nucleation rate of type-i bubbles. In these equations we are neglecting 
"secondary" bubble nucleation which may occur within a comoving distance from any 
given "primary" nucleation event, before the primary bubble reaches its asymptotic cooving 
size H~^. Secondary nucleations will affect the co-moving volume distribution very little, 
especially if the nucleation rate per unit volume is much smaller than H'^. This is a typical 
situation since the nucleation rates are usually exponentially suppressed. The solution of 
these equations with the initial condition (|T^ is 



1 — exp 



47rA 



(20) 



where A = Ai + A2. The fraction of markers that end up in bubbles of type i is given by 

Pi = fi{t ^ 00) = A,/A. (21) 
Now, substituting (|T6|) and (^) into Eq.® for the probabilities we obtain 



-Pi _ Ai / Z^i ^ 

P2 A2 \Z^2' 



(22) 



This agrees with one's intuitive expectation that the probability should be proportional to 
the nucleation rate and to the volume expansion factor inside the bubbles. 



V. DISCUSSION 

In this paper we have suggested a cutoff procedure which allows one to assign probabilities 
to different types of thermalized regions in an eternally inflating universe. The probabilities 
are calculated with the aid of "markers" - imaginary pointlike objects which are assumed 
to be created at a constant rate in the inflating regions where the infiaton field is at the top 
of its potential. The probability for regions of type i is then 

P^ cx p^N,, (23) 

where pi is the fraction of markers that end up in type-i regions and Ni is the number of 
galaxies formed in a comoving sphere of radius equal to the average separation between the 
markers. 

In contrast to some earlier proposals, the new prescription is manifestly gauge- invariant. 
It also gives the expected results in cases where we have well-motivated intuitive expectations 
for the probabilities. 

The method of calculating probabilities presented in this paper has some similarities 



to the so-called e-prescription which was proposed in Ref. [0]. Starting with a comoving 
volume with near the top of the potential, the numbers of galaxies are calculated in this 
prescription by imposing cutoffs at different times in different types of regions. The cutoff 
in type-i regions is chosen at the time td when all but a small fraction e of the comoving 
volume destined to thermalize in this type of regions has thermalized. The value of e is the 
same for all types of regions, but the cutoff times are different. The limit e — >^ is taken 
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after calculating the probabilities. The e-prescription was applied in to calculate the 
probabilities in open inflation and gave the same result (p2|) that we obtained here. 

To see the connection between this prescription and the method of the present paper, 
imagine that the initial comoving region contains a large number of markers. If Pi is the frac- 
tion of markers that are going to end up in type-? regions, then the cutoffs in e-prescription 
are imposed at the times when the numbers of markers N'mi in the two types of regions are 
related by N'mi/f^m2 = Pi/P2- The number of galaxies in each type of region is Afgi = AfmiNi, 
where Ni is the number of galaxies per one marker. Hence, 

which has the same form as (||). 

One difference between the two methods is that markers are continuously produced in 
our new approach, while in e-prescription new markers are not produced even in regions 
when the infianton field fluctuates back to the top of its potential. Another difference is 
that e -prescription uses constant-t cutoffs, while the new approach uses spherical cutoffs. 
Because of its reliance on a constant-t cutoff, the results obtained using e-prescription are 
generally gauge- dependent, whereas the new method is gauge independent. 

The most straightforward way to implement the new prescription is through a numerical 
simulation of an eternally inflating spacetime. This method, however, suffers from severe 
computational limitations |jTl||. Alternatively, one can use an approximate analytic method 



based on the Fokker-Planck (FP) formalism of stochastic inflation []7|,pT|,^ . This method 



works very well for the calculation of Pi . However, the results for the density of markers 
obtained using this approach are generally not gauge invariant. 

The FP formalism can be used to calculate the physical volume that thermalizes 

prior to some time t and the number of markers N^i{t) contained in that volume for different 
types of regions. The average distance between the markers can be expressed as 



di{t) 



One might expect that in the limit t — > 00, di{t) should approach the average separation 
between markers on the thermalization surfaces (calculated in large spherical volumes). 
However, this is not generally the case. The thermalized volume Vn{t) in an eternally 
inflating universe is dominated by the newly thermalized regions (shortly before the time 
t), and the density of markers in these regions is generally correlated with the choice of the 
time variable t. 

The density of markers on the surfaces can vary greatly from one area to another. 
In some places the field takes a more or less direct route from the top of the potential 
to resulting in a relatively high density of markers, while in others it takes a long time 
fluctuating up and down the potential, so that the markers are greatly diluted. By including 
only regions that thermalize prior to time t, we "reward" regions that thermalize faster and 
therefore introduce a bias favoring higher densities of markers. This qualitative tendency 
is present for most choices of the time variable, but quantitatively the bias will not be the 
same. Hence, one should not be surprised that di calculated from Eq. (|25|) would depend on 
the choice of gauge. 
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Despite this gauge dependence, the FP method may give approximately vahd results for 
some class of models. It has been argued in Refs. |19,23|| that in the case of e-prescription 
the gauge dependence is rather weak for a wide class of potentials. One can expect the 
situation to be similar for our new method when a constant t cutoff is used. However, more 
work is needed to determine what additional requirements the potential should satisfy for 
this approximate gauge independence to apply. 

We started in Section II of this paper by formulating a set of requirements which we 
believe any method for calculating the probabilities should satisfy. The specific prescription 
we introduced in Section III can be regarded as an "existence proof", demonstrating that 
a prescription satisfying all the criteria can indeed be constructed. It is quite possible, 
however, that our prescription is not unique and that more attractive and better motivated 
methods can be developed. With this in mind, we conclude by indicating some possible 
shortcomings and limitations of our method. 

Admittedly, our prescription includes an element of arbitrariness when we assume that 
markers are produced only in regions where the inflaton is at the top of its potential. 
This "5-function" source should not be understood literally. The semiclassical picture of 
eternal inflation involves smearing over spacetime scales ~ and over scalar field intervals 
~ H/27T. Hence, marker formation at the top of the potential is equivalent to marker 
formation within an interval A0 ~ H/2tc from the top of the potential. However, an attempt 
to extend marker formation to a wider interval encounters some ambiguities. Formally, there 
is no problem in calculating the density Ui of markers in thermalized regions even if these 
are produced at some given rate R{(l))d(f) per unit proper time and volume in regios where 
the field is in the range d(j). However, there is an obvious arbitrariness in what should be 
chosen as our "smearing function" R{(j)). Also, if markers are produced at different rates at 
different values of H{(j)), the calculation of the fraction of markers Pi that end up in region 
i becomes somewhat ill posed. Markers formed at different values of (p will have a different 
probability pi of ending up in a type-i region and it is not clear how to weigh the different 
contributions. One could calculate Pi{(f>) and then average over with some weight, but it is 
not clear how the weight function is to be determined. It thus appears that confining marker 
formation to the top of the potential has some advantages and may be not as arbitrary as 
it first seems. 

It is not clear whether or not our method will give reasonable results when applied to 
the most general type of inflationary scenarios, when the inflaton potential has several local 
maxima. Additional complications arise in models where the minimum of the potential 
has a no n- vanishing cosmological constant. In such models, regions of true vacuum in the 
post-inflationary universe may fluctuate back to the quantum diffusion range of the inflaton 
potential and the spacetime structure is more complicated than the one represented in Fig. 2. 

Also, our method is not applicable when the potential is unbounded from above, in which 
case eternal inflation runs into the Planck boundary. This is not particularly worrisome, since 
after all the inflaton may be a modulus with a compact range, and the inflaton potential 
may well be bounded from above at a scale much lower than the Planck scale. 
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